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Abstract 

Background: Rapid detection and therapeutic intervention for infectious and emerging diseases is a major scientific 
goal in biodefense and public health. Toward this end, cytokine profiles in human blood were investigated using a 
human whole blood ex vivo exposure model, called WEEM. 

Results: Samples of whole blood from healthy volunteers were incubated with seven pathogens including Yersinia 
pseudotuberculosis, Yersinia enterocolitica, Bacillus anthracis, and multiple strains of Yersinia pestis, and multiplexed 
protein expression profiling was conducted on supernatants of these cultures with an antibody array to detect 30 
cytokines simultaneously. Levels of 8 cytokines, IL-1 a, IL-1(3, IL-6, IL-8, IL-10, IP-10, MCP-1 and TNFa, were significantly 
up-regulated in plasma after bacterial exposures of 4 hours. Statistical clustering was applied to group the 
pathogens based on the host response protein expression profiles. The nearest phylogenetic neighbors clustered 
more closely than the more distant pathogens, and all seven pathogens were clearly differentiated from the 
unexposed control. In addition, the Y. pestis and Yersinia near neighbors were differentiated from the B. anthracis 
strains. 

Conclusions: Cluster analysis, based on host response cytokine profiles, indicates that distinct patterns of 
immunomodulatory proteins are induced by the different pathogen exposures and these patterns may enable 
further development into biomarkers for diagnosing pathogen exposure. 



Background 

Yersinia pestis and Bacillus anthracis are two pathogens 
of significant concern to public health from a biodefense 
perspective [1,2]. Y pestis, the causative agent of plague, is 
a Gram-negative, highly communicable coccobacillus that 
has been responsible for three historic pandemics with 
high mortality rates [3-5]. The microorganism possesses a 
Type III secretion mechanism common to several human, 
animal and plant pathogens, whereby a series of pathogen- 
specific structural proteins form a syringe-like structure 
capable of injecting virulence factors into the mammalian 
host cell. These virulence factors then facilitate pathogen 
use of host nutrients and thwart the host immune response, 
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ultimately causing cell and host death [6,7]. Naturally occur- 
ring plague can be transmitted from infected fleas and 
rodents to humans, and although the pathogen can be pha- 
gocytosed, it can also resist destruction by manipulating the 
host defense mechanism(s), potentially through antigenic 
mimicry [8]. Y pestis then multiplies rapidly leading to ne- 
crosis of lymph nodes, a condition known as bubonic 
plague, which can result in death if untreated [2]. In some 
cases the infection can spread through the blood stream 
resulting in systemic plague (septicemia) or to the lungs 
resulting in the highly contagious and deadly form of the 
disease known as pneumonic plague. There are currently 
no rapid, widely available diagnostic tests for plague, and 
the most common treatment is streptomycin [2,3], an anti- 
biotic with adverse effects. 

Two other species from the genus Yersinia are also 
human pathogens: Y pseudotuberculosis and Y 
enterocolitica [9,10]. Despite their high degree of se- 
quence similarity to Y pestis, these two near neighbors 
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of Y. pestis manifest in very different symptoms, ran- 
ging from abdominal pain to septicemia in humans, 
usually caused by infection through contaminated 
food. Infections caused by Y. pseudotuberculosis or Y, 
enterocolitica can be effectively treated with antibiotics 
and in most cases are self-limiting. Notably, Y pestis is 
reported to have evolved from Y pseudotuberculosis 
within the past 10,000 years [11]. 

B. anthracis is a Gram-positive, rod-shaped spore- 
forming bacterial pathogen and the causative agent of 
anthrax [12,13]. Human, livestock, and wildlife mortal- 
ities attributable to anthrax occur in numerous regions 
of the world, although the majority of cases are found 
in less industrialized nations [14]. Three forms of the 
disease have been described: cutaneous, intestinal and 
inhalational. While cutaneous and intestinal forms may 
be less severe, inhalational anthrax is often fatal without 
prompt antibiotic treatment [13]. The primary mechan- 
isms of virulence employed by B. anthracis are asso- 
ciated with two virulence plasmids designated pXOl 
and pX02 [15]. The net effect of these plasmids is vir- 
tually unhindered proliferation of B. anthracis within 
the host, hemorrhaging, cardio-pulmonary collapse, and 
death. 

The regulation of production of host cytokines by both 
Yersinia and B. anthracis has been described previously. 
Pickering A. K. et. al. measured cytokine levels in human 
dendritic cell supernatant and in mouse peritoneal 
macrophages exposed to B. anthracis spores [16]. They 
observed significant increase in TNF-a, IL-6, IL-lp, IL-8, 
and IL-12 in human dendritic cell supernatants by 5 
hours post-exposure. High levels of IL-6, and TNF-a 
were observed in the supernatant from B. anthracis 
infected mouse peritoneal macrophages [16]. In a mouse 
model, 6 cytokines, namely IL-12p70, TNF, IFN-y, MCP-1, 
IL-10, and IL-6, were increased significantly in mouse lung 
at 48 hours of Y. pestis infection [17]. In previous work 
comparing exposures to different bacterial pathogens, 
distinct patterns of cytokine expression levels were found 
that could discriminate the particular host response [18], 
including while using pathogen-specific LPS in whole 
blood [19]. 

The hypothesis for the present study is that ex- 
posure to diverse bacterial pathogen strains would 
result in distinct cytokine profiles in the host, with 
strains from the same species exhibiting more simi- 
lar profiles than strains from phylogenetically dis- 
tant species. A multiplex cytokine protein chip was 
used, and a multivariate approach was taken that 
combined expression data on multiple cytokines. 
Multivariate clustering techniques were used to es- 
tablish cytokine expression profiles after ex vivo ex- 
posure of whole blood to seven pathogens. 



Methods 

Bacterial strains and culture conditions 

The bacterial strains used in this study include: B. 
anthracis Ames (virulent), B. anthracis Sterne (vaccine 
strain), Y. pestis KIM5 D27 (attenuated, pgm-). Y. pestis 
India/P (attenuated, pgm-), and Y. pestis NYC (virulent), 
Y. pseudotuberculosis serotype 1 PB1, and Y. enterocoli- 
tica WA serovar 0:8. Bacteria were grown on tryptose 
blood agar slants at 26°C for 1-2 days and subsequently 
collected using 2 ml of 0.033M potassium-phosphate, 
pH 7.0;.bacterial densities were measured at OD 620 
(1 OD 620 = L2 x 10 9 colony forming units/ml). 

Whole blood ex vivo exposure model (WEEM) 

Human blood was collected from a healthy donor by 
venipuncture using CPT Vacutainer tubes (Becton Dick- 
inson) containing citrate. Informed consent was obtained 
and our blood collection protocol was approved by the 
LLNL IRB committee. Separate CPT tubes were used for 
the unexposed control and 7 different bacterial expo- 
sures (B. anthracis Ames, B. anthracis Sterne, Y. pestis 
NYC, Y. pestis India/P, Y. pestis KIM5 D27, Y. pseudotu- 
berculosis, and Y. enterocolitica). Bacteria were added to 
blood within 15 minutes of collection at a multiplicity of 
infection ratio of 5:1. This ratio was determined against 
white blood cells in whole blood:~7 x 10 6 cells/ml. Each 
whole blood sample was incubated with bacteria for 4 
hours at 37°C in 5% C0 2 Following incubation, plasma 
was collected by centrifugation at 2000 x g for 10 min at 
4°C. The control plasma was obtained in the same way 
and treated with 0.033 M potassium-phosphate as a 
mock exposure. These plasma samples were used for 
cytokine measurements. 

Cytokine immunoassays with protein arrays 

The measurements of cytokines were performed using 
Zyomyx Protein Profiling Biochips (Hayward, CA). These 
protein arrays allow the simultaneous quantification of 
30 biologically relevant cytokines, as determined by Zyo- 
myx, Inc: IL-la, IL-lp, IL-2, IL-3, IL-4, IL-5, IL-6, IL-7, 
IL-8, IL-10, IL-12(p40), IL-12(p40/p70), IL-12(p70), IL- 
13, IL-15, TNFa, TNF (3, Eotaxin, MCP-1, MCP-3, 
TRAIL, CD95(sFas), MIG, sICAM-1, IP-10, CD23, TGF- 
p, GM-CSF, GCSF, IFN-y. Each cytokine assay was opti- 
mized for the Zyomyx Protein Profiling Biochip based on 
many factors including the availability of antibodies and 
the sensitivity and specificity of antibody-cytokine inter- 
actions. Each protein array chip is designed with 6 inde- 
pendent microfluidic channels that allow up to 6 
samples to be loaded into isolated regions of an array. 
Antibodies specific for 30 analytes were arrayed in each 
channel, and each antibody was arrayed in redundancy 
on 5 pillars within the channel. Accordingly, a cytokine 
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measurement represents the average of 5 measurements. 
All immunoassay steps, including sample loading, wash- 
ing, and detection, were performed with a fully auto- 
mated biochip processing station (Zyomyx Assay 1200 
workstation). Eight protein array chips were used in 
these experiments. Two chips were used for generating 
calibration curves with a calibration standard kit con- 
taining 30 analytes (Zyomyx, Inc.). Sample (40 ul) was 
injected into each channel of the protein array chips. 
Standard solutions were applied to two channels of each 
chip for chip-to-chip normalization. Triplicates of con- 
trol and pathogen-exposed plasmas were applied ran- 
domly to four channels of 6 protein array chips. Protein 
arrays were scanned at 532 nm with Zyomyx Scanner 
100 after immunoassays. Zyomyx Data Reduction soft- 
ware was used for normalization, calculation of calibra- 
tion curves. Dixons test was used to remove outliers, 
and the median feature intensity was background sub- 
tracted. Concentrations of cytokines in plasma samples 
were determined by a four parameter logistic model. 

Cluster analysis of cytokine data 

Multiple hierarchical clustering methods were used to 
group the pathogen exposures based on the multivariate 
cytokine expression profiles induced in a host infection 
model system. First, hierarchical agglomerative clustering 
[20] was applied to group the control and the seven 
pathogen-exposed samples based on their cytokine con- 
centration profiles. Each of the eight samples was charac- 
terized by the multivariate vector of its average loglO 
concentration over the cytokines. The samples were then 
clustered based on the following distance measures be- 
tween the samples and between the clusters. 

Distance between two samples was defined using two 
distance metrics: 

• Euclidean distance 

• Correlation distance: (1 - Spearman correlation 
coefficient between the samples) 

Distance between two clusters was defined using three 
methods: 

• Complete linkage (furthest neighbor): the largest 
distance between members of the clusters 

• Single linkage (nearest neighbor): the smallest 
distance between members of the clusters 

• Average linkage (group average): the average 
distance between members of the clusters 

Given a pair of distance metrics between samples and 
clusters, the algorithm was initialized with the eight sam- 
ples forming eight different clusters and then processed 



iteratively by joining the two most similar clusters. The 
tree was built starting from the individual samples, using 
an agglomerative (bottom up) approach. The resulting 
hierarchy of clusters was displayed as a dendrogram. 
These traditional clustering methods provide a quick, ex- 
ploratory overview of the data. However, these methods do 
not estimate the optimal number of clusters in the data; ra- 
ther, the clustering is performed exhaustively from the low- 
est possible level of the hierarchy where each sample forms 
its own cluster, to the highest level where all samples are 
grouped into one cluster. 

In addition to the traditional hierarchical agglomera- 
tive clustering method, the hierarchical ordered parti- 
tioning and collapsing hybrid (HOPACH) algorithm was 
also applied to the cytokine measurements [21]. In con- 
trast with the previous approaches where the tree was 
built starting from the individual samples as the leaf 
nodes, HOPACH used a hybrid divisive-agglomerative 
approach: it started from the root cluster containing all 
the samples (divisive, top down approach), then divided 
the root down to leaf nodes, with an extra collapsing (ag- 
glomerative) step after each iteration that combined simi- 
lar clusters. Based on the correlation distance between 
samples, HOPACH determined the split that minimized a 
measure of cluster homogeneity called the median split sil- 
houette. While computationally more expensive than the 
previous methods, HOPACH was expected to perform bet- 
ter because of its dynamic approach to update and poten- 
tially revise the clusters at every step of the iteration. 
Furthermore, HOPACH also estimated the optimal num- 
ber of clusters from the data, and thus offered another ad- 
vantage over the previous methods. 

Computations were performed in the R computing en- 
vironment (http://www.r-project.org/) and the HOPACH 
package [21]. 

Results 

Cytokine levels were examined using an ex vivo model, 
termed WEEM for whole blood ex vivo exposure model. 
Individual samples of anti-coagulated human blood were 
incubated with B. anthracis Ames, B. anthracis Sterne, Y 
pestis KIM5 D27, Y. pestis NYC, Y. pestis India/P, Y 
enter vcolitica, and Y. pseudotuberculosis. After 4 hours 
exposure, blood cells were removed by low-speed centri- 
fugation and concentrations of 30 cytokines in the 
plasma were measured with protein arrays. Concentra- 
tions of fourteen cytokines, GCSF, IFNy, GM-CSF, IL-7, 
IL-12(p70), IL-12(p40/p70), IL-13, IL-2, IL-3, IL-4, IL-5, 
MCP-3, TGFp, and TNFp were below the limit of detec- 
tion in this study. The following 16 cytokines were 
detected: Eotaxin, IL-10, IL-12(p40), IL-15, IL-la, IL-1(3, 
IL-6, IL-8, IP-10, MCP-1, MIG, TNFa, TRAIL, sCD23, 
sCD95, and sICAM-1 (Figure 1). To determine if there 
were significant differences among the levels of cytokines 
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Figure 1 Scatter plots of 16 cytokine concentrations detected in human blood following ex vivo bacterial exposures. Cytokine 
concentrations were displayed on a logarithmic scale. The cytokines shown here were detected out of the 30 cytokines in the arrays. The 8 
cytokines that were found to be statistically differentially expressed among these samples are highlighted with rectangular boxes. Each mark 
delineates the average of triplicate exposure samples. Each exposure sample is loaded onto a protein array chip that contains 5 independent 
measurements per cytokine meaning that fifteen measurements are used to obtain these data. 



in the control and pathogen exposed plasma samples, F- 
tests were performed. For thirteen of these 16 cytokines, 
all three replicates were detected and these cytokines 
were subjected to F-tests. Statistical analysis indicated 
that 8 cytokines (IL-la, IL-lp, IL-6, IL-8, IL-10, IP-10, 
MCP-1, and TNFa) had differentially elevated expression 
profiles following different bacterial exposures. Figure 2 
shows the concentrations (pg/ml) of these cytokines in 
the control and bacteria exposed plasma samples. The F- 
tests revealed that the other five cytokines containing 
complete datasets, TRAIL, sCD23, sCD95, MIG, and 
sICAM-1, had no significant difference between bacterial 
exposures and the mock-exposed control. Moreover, 
there was a great variation in absolute concentrations be- 
tween cytokines. For example, the concentrations of 
TNFa, sCD23, and sICAM-1 were as high as 1 x 10 4 
-10 5 pg/ml, whereas IL-10 was much lower, about 16 pg/ 
ml. 

Marked differences in induced cytokine patterns be- 
tween B. anthracis and Yersinia exposures were found. 
Also, the levels of induction of these cytokines differed 
among the different bacteria. For example, Yersinia spe- 
cies induced much higher cytokine response than B. 
anthracis for IL-la, IL-lp, IL-6, and TNF-a (Figure 2). 
The two strains of B. anthracis bacteria induced different 
levels of IL-1(3 and TNF-a (Figure 2), including 2.2 times 
higher concentration of IL-lp and 1.6 times higher con- 
centration of TNF-a for the Sterne strain than the Ames 
strain of B. anthracis. These differences were statistically 
significant (pairwise £-test p value = 0.0039 for IL-lp and 
0.022 for TNF-a). To discriminate Y. pestis exposure 



from near neighbors, IL-10 levels can be used, showing 
cytokine concentrations following Y. enterocolitica ex- 
posure and Y pseudotuberculosis exposure that are on 
average 5-fold higher and 2-fold higher, respectively, 
than after Y. pestis exposure (Figure 2). IL-10 differential 
expression was specific to the Yersinia spp. because ex- 
posure to B. anthracis strains showed comparable IL-10 
levels to that in unexposed control. 

The HOPACH algorithm estimated the number of 
clusters as five, and grouped the samples based on their 
host cytokine expression profiles as follows: 1) Y. pestis 
(KIM5 D27, India/P, and NYC), 2) Y pseudotuberculosis, 
3) Y enterocolitica, 4) B. anthracis (Ames and Sterne), 
and 5) Control (Figure 3). The closer the pathogen- 
exposed samples are within the tree on the left, the more 
similar they are. Height of the branches indicates the dis- 
tance between the successive nodes in the clustering. 
The method separated the B. anthracis and Yersinia 
infected blood samples. In addition, the cytokine profile 
of the mock-exposed control was more similar to the 
pattern produced by B. anthracis exposure than to the 
profile elicited by Yersinia. 

Results of the hierarchical clustering when using the 
Euclidean distance between samples depended on the 
distance metric between clusters. The three methods for 
determining the distance between clusters (complete 
linkage, single linkage, and average linkage, see Materials 
and Methods) all established three major clusters: 1) Y 
pestis and near neighbors, 2) B. anthracis, and c) Con- 
trol. Differences between the results occurred when the 
Yersinia cluster was further divided. The average linkage 
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Figure 2 Concentrations of 8 cytokines in human whole blood after ex vivo exposure to pathogens. The control was a mock-exposed 
sample. Cytokine concentrations were determined using protein arrays. The bars represent the average of three replicate samples that each 
contain 5 replicate features per cytokine assay and the lines represent the standard deviation among the three replicates. 



method, consistent with Figure 3, formed a subgroup of 
the three Y. pestis strains, then grouped them first with 
Y. pseudotuberculosis followed by Y. enter ocolitica. 
Complete and single linkage methods, however, first 
grouped the attenuated virulent strain of Y pestis (India/ P) 
with the more virulent strain (NYC), both clinical isolates 
from human plague cases, and then clustered them with 
Y pseudotuberculosis, followed by the attenuated Y pestis 
(KIM5 D27), and lastly with Y enter ocolitica. This is inter- 
esting from an evolutionary perspective because it has been 
proposed that Y pestis evolved from Y pseudotuberculosis 
within the last 10,000 years, and thus these two pathogens 
are more closely related [11]. 

When using hierarchical clustering with the correlation 
distance between the samples, the final clusters were in- 
dependent of the distance metric between clusters, and 
agreed with the tree structure in Figure 3. The complete, 



single, and average linkage methods all resulted in the 
following major clusters: 1) Yersinia, 2) B. anthracis, and 
3). Control. Within the Yersinia cluster, Y. pestis (NYC) 
was closest to Y pestis (India/P), followed by Y pestis 
(KIM5 D27), Y pseudotuberculosis, and Y enter ocolitica. 

Discussion 

The HOPACH clustering method (Figure 3) produced 
five distinctly separated clusters: 1) Y pestis (KIM5 D27, 
India/P, and NYC), 2) Y. pseudotuberculosis, 3) Y. enter o- 
colitica, 4) B. anthracis (Ames and Sterne), and 5) Con- 
trol. This result is consistent with the findings using the 
correlation distance and the Euclidean distance with 
average linkage. In addition, HOPACH estimated the op- 
timal number of clusters as five. That is, the Yersinia 
subcluster is best if it is divided into the three clusters 
specified by 1) through 5) above. Y. enterocolitica forms 
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Figure 3 Clustering result with HOPACH using the average linkage distance between clusters is shown. The eight pathogen-exposed 
samples are clustered according to the dendrogram on the left and cluster into five groups, 1) Y. pestis (KIM5, NYC, and India), 2) Y. 
pseudotuberculosis, 3) Y. enterocolitico, 4) B. anthracis (Ames and Sterne), and 5) Control. Sixteen cytokines (Eotaxin, IL-10, IL-12(p40), IL-15, IL-1a, IL- 
1(3, IL-6, IL-8, IP-10, MCP-1, MIG, TNFa, TRAIL, sCD23, sCD95, and slCAM-1) are also reordered based on their correlations according to the 
dendrogram on the top. Clusters go from root at top to leaf node for each cytokine. Clusters in between are based on their agglomerative . The 
branch shows the similarity, the short the branch, the more similar. In addition, the eight rightmost proteins form a cluster that may involve 
inflammation-related cascades initiated by an innate immune response to these pathogen. Colors represent units of log 10 [pg/ml], in ten equally 
spaced intervals increasing from white to dark red. A key showing the specific log 10 values for each interval is shown in the figure. 



its own cluster, and so does Y. pseudotuberculosis. Y. pes- 
tis (KIM5 D27), Y pestis (India/P), and Y pestis (NYC) 
are grouped into one cluster. Further subdivisions lead 
to an overall clustering with inferior quality. 

In addition to clustering the cytokine expression pro- 
files across bacterial treatments, Figure 3 also groups the 
cytokines themselves and clusters the proteins based on 
their similarities across the pathogen exposures and reor- 
ders them accordingly. Interestingly, the three pro-in- 
flammatory cytokines IL-lp, TNFa, and IL-6 clustered 
closely, and so did the three chemokines MCP-1, IP-10, 
and IL-8. Although these 6 cytokines do not cluster as a 
single group, they do cluster at a branch further away 
from the leaf node, which includes IL-10 and sCD95, to 
make a larger group of 8 proteins. Several of these pro- 
teins are involved in inflammatory conditions, such as 
IL-lbeta, TNFa, IL-6, [22] and have been shown to be 
upregulated in cell culture and animal model specifically 
exposed to biothreat agents [23]. Increased expression of 
IL-6 and TNFa clustered together in a study involving 
mouse splenic CDllb + cells following sub-lethal Y 
enterocolitica infection [24]. In addition, several of the 
cytokines in this cluster, namely TNF-alpha, ILl-beta, 
IL-10, and MCP-1 are expressed higher in exposed whole 
blood as compared to control in this study and in whole 
blood exposure to LPS from several other gram negative 
bacterial pathogens [19]. In addition to expression differ- 
ences, the absence of detected cytokine expression can 
also be helpful in discriminating pathogen exposure. 



The multiplex detection of 30 cytokines in this study 
revealed the early phase cytokine expression profiles in 
human plasma following exposures to B. anthracis 
(Ames and Sterne), Y. pestis (KIM5 D27, NYC and 
India/P), Y. pseudotuberculosis, and Y. enterocolitica. 
The expression levels of 8 cytokines, IL-la, IL-lp, IL-6, 
IL-8, IL-10, IP-10, MCP-1, and TNFa were significantly 
different from that of unexposed control (Figure 2). Al- 
though the focus of our work was to show that cytokine 
expression profiling can discriminate between different 
pathogen exposures in a human whole blood ex vivo 
model, these results also represent an initial attempt to 
characterize the full cytokine response to each individual 
pathogen. Our preliminary study using a single exposure 
protocol at a single time post-exposure will need to be 
supplemented with more thorough investigation in order 
to determine the usefulness of using cytokine levels for 
diagnosing pathogen exposure. However, the single time 
point chosen, 4 hours, is sufficient to detect proteomic 
changes and has been used in previous studies examin- 
ing cytokine levels [25-27]. This time point represents a 
start towards a more complete temporal study\, as has 
been done with gene expression patterns for two of the 
pathogens studied here [25,27]. In addition, studies that 
provide expression patterns for a single cytokine using 
multiple time points will also be needed to make the 
results of this paper clinically useful, such as has been 
done by, Cooper and coworkers, who examined IL- 
12p40 and IL-12p70 levels following different growth 
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conditions and exposure levels for a time course of Y 
pestis exposed dendritic cells [28]. The results of the 
current work shows a similar expression pattern trend 
to this previous work, in which, Y pestis induces IL- 
12p40 and at a substantially higher level than IL-12p70. 

Our results showed that the expression levels of 3 
chemokines, IL-8, MCP-1 and IP- 10, were induced by 
both Yersinia and B. anthracis exposures. No significant 
differences were found for these cytokines between 
Yersinia and B. anthracis exposures. IL-8, MCP-1 and 
IP- 10 are chemokines that enable the migration of 
leukocytes from the blood to the site of inflammation. 
IL-8 is a key chemokine regulating neutrophil recruit- 
ment [29]. The essential involvement of IL-8 in acute 
inflammation was demonstrated by neutralizing IL-8 
with its antibody. When highly specific antibody against 
IL-8 was administered in acute inflammatory reactions 
induced by several stimuli including lipopolysaccharide, 
neutrophil infiltration was blocked [30]. MCP-1 is 
known for its ability to act as potent chemoattractant 
and activator of monocytes/macrophages as well as NK 
cells but not neutrophils [31,32] . IP-10 has no chemo- 
tactic activity for neutrophils but attracts monocytes, 
NK, and T cells to the site of infection and regulates T 
cell maturation [33,34]. It was reported previously that 
elevated IL-8 and MCP-1 were secreted by human epi- 
thelial cells after Y. enter ocolitica infection, but not IP- 
10 [35,36]. Human dendritic cells, infected with B. 
anthracis spores, secreted high level of IL-8 at 7.5 hours 
[16]. In our study, the fold increase of IL-8 was much 
greater than MCP-1 and IP-10 (Figure 2). For example, 
the induction of IL-8 by Ames strain of B. anthracis 
was 41 fold, while MCP-1 was 2 fold and IP-10 was 2.5 
fold (Figure 2). This result may indicate that IL-8 is a 
dominant chemokine in early response (4 hours expos- 
ure in our study) and neutrophils are the major player 
in early inflammatory response. 

Here we compared cytokines induced by B. anthracis 
and Yersinia exposures. Overall, Yersinia exposure 
induced higher levels of IL-la, IL-1(3, IL-6, IL-10 and 
TNFa than B. anthracis exposure, suggesting these cyto- 
kines could be used to develop an assay for discriminat- 
ing Yersinia spp. from B. anthracis exposures. The 
vaccine strain (Sterne) of B. anthracis induced higher 
levels of IL-lp and TNFa than the virulent strain (Ames) 
(Figure 2), suggesting these cytokines can contribute to a 
biomarker panel to discriminate if a particular isolate of 
B. anthracis is virulent. There was also a difference in in- 
duction of IL-10 between Y. pestis and near neighbors 
(Figure 2), suggesting this cytokine is a candidate bio- 
marker for discriminating the virulence of Yersinia spe- 
cies. These data regarding IL-10 expression following 
Yersinia spp. exposure are in agreement with published 
literature that shows Y enter ocolitica and Y pestis can 



elicit statistically different levels of IL-10 expression [37]. 
Differences in IL-10 induction may be due to differences 
in the lcrV protein among Yersinia spp. [38]. The differ- 
ent cytokine profiles induced by B. anthracis and Yersi- 
nia here may be partially due to different surface 
antigens on the outermost part of these pathogens and 
the manner in which these bacteria were grown. Lipo- 
polysaccharide (LPS), the main constituent of the outer 
membrane of Gram-negative bacteria, and peptidoglycan 
(PGN), the major cell wall component of Gram-positive 
bacteria, have been reported to elicit markedly different 
immune responses [39]. However, virulence factors, such 
as B. anthracis lethal toxin and Yersinia virulence anti- 
gen, LcrV, may also play important roles in differential 
cytokine induction. This view is supported by numerous 
reports that B. anthracis toxin and virulence factors of 
Yersinia bacteria (Yops, invasin, LcrV) modulate host 
cytokine responses [40-51]. 

While the various clustering methods resulted in 
slightly different final hierarchies, all were consistent in 
separating the unexposed control from the samples 
exposed to B. anthracis or to the Y pestis and near 
neighbors. Agreement on this level among the various 
clustering procedures lends more confidence to the over- 
all results. On a more detailed level, the methods 
grouped slightly differently the samples exposed to the Y. 
pestis and near neighbors, which indicates that these 
samples cannot be unequivocally separated based on the 
current data and additional biomarkers or a larger sam- 
ple set would be needed. The most advanced HOPACH 
method estimated the optimal number of clusters in the 
data as five, corresponding to the unexposed control, 
and the four species: B. anthracis, Y. pseudotuberculosis, 
Y enter ocolitica, and Y pestis (avirulent and virulent) 
(Figure 3). 

Information gained from the targeted protein array 
data for host response complements genomic [52-56], 
and other proteomic studies [57-60] of host-pathogen 
interactions. The success of the WEEM and computa- 
tional method to distinguish pathogen exposure, based 
on host response in this initial study, is encouraging and 
suggests a number of possibilities for future studies to 
refine the findings. Comparative analysis, such as the 
current work, can potentially reveal the critical patho- 
genic mechanism(s) and host innate immune responses 
during infection as was previously shown for Y. pestis 
and Y. pseudotuberculosis [61]. Opportunities include 
using statistical hypothesis tests based on analysis of 
variance to assess the significance of the observed differ- 
ences among the host-pathogen cytokine concentration 
profiles, as well as performing follow-up studies to focus 
more on the Y pestis and near neighbor cluster. In 
addition, the methods can be extended to investigate 
host responses to diverse pathogens in multiple host 
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model systems to cross validate the significance of the 
biomarkers to distinguish pathogen exposures. 

Conclusion 

Results from this study suggest that cytokine arrays 
coupled with statistical clustering methods can distin- 
guish exposures to pathogens, including multiple strains 
of Y. pestis, Y. pseudotuberculosis, Y. enterocolitica, and 
B. anthracis. These methods differentiate both near 
neighbors and distant evolutionary microbes based on 
host response data. The distinct cytokine profiles also 
provide insight into both the host response and virulence 
mechanisms of diverse pathogens. In summary, 
characterization of host responses based on cytokine 
profiles has translational application, potentially provid- 
ing the identification of infectious diseases and leading 
toward the ultimate goal of presymptomatic detection 
via sentinel surveillance of pathogen exposure and ap- 
propriate treatment. 
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